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1.0  INTRODOCTION 

Large-amplitude,  low-frequency  pressure  oscillations  have  been  recognized 

1  -9 

as  a  major  limiting  factor  in  the  performance  of  ramjet  propulsion  systems. 

The  oscillations  arise  from  combustion  instabilities  occurring  in  the 
combustor,  and  from  the  inherent  instability  of  the  transonic  diffuser  in  the 
inlet  system.  Coupling  between  these  two  systems  can  lead  to  unacceptably 
high  pressure  levels  which  degrade  system  performance,  unstart  the  inlet,  and 
cause  structural  failure. 

McDonnell  Douglas  Research  Laboratories  (MDRL)  has  been  examining  inlet 
system  unsteadiness  through  studies  of  the  self-excited  and  externally  excited 
flowfields  in  transonic  diffusers.  The  flowfields  fall  into  three  basic 

classes,  as  depicted  in  Fig.  1.  The  type  of  unsteadiness  observed  depends 
strongly  on  the  type  of  mean  flowfield  which  exists  in  the  diffuser. 

Moderate- area-ratio  diffusers  with  relatively  weak  shocks  (Fig.  la)  contain 
fully  attached  flow  with  thin  boundary  layers.  A  well-defined  core  flow 
exists  to  the  end  of  the  diffuser.  For  weak  shocks  in  large-area-ratio 
diffusers  (Fig.  1b),  a  pressure-gradient-induced  separation  may  occur. 

However,  the  boundary  layers  remain  relatively  thin.  For  all  strong-shock 
flows,  independent  of  area  ratio  (Fig.  1c),  the  flow  exhibits  a  shock-induced 
separation.  The  boundary  layers  grow  rapidly  and  merge  within  the  diverging 
section,  terminating  the  core  flow. 

Three  fundamental  mechanisms  have  been  proposed  as  contributors  to  the 
1 3 

observed  unsteadiness  (Fig.  2):  upstream  and  downstream- traveling  acoustic 
waves,  and  a  disturbance  which  is  convected  through  a  thick,  separated 
boundary  layer,  such  as  depicted  in  Fig.  1c.  The  convective  wave  is  generated 
by  the  recoil  of  the  shock  after  impingement  of  an  upstream-traveling  acoustic 
wave . 

A  common  method  used  to  examine  unsteady  flowfields  in  a  controlled  manner 
is  to  subject  the  flow  to  a  known,  periodic  excitation.  This  method  has 

proven  useful  in  the  study  of  coherent  structures  in  free  shear  layers  and 

19-21  22-23 

jets,  as  well  as  in  air-breathing  propulsion  systems.  This 

technique  allows  the  overall  response  of  the  system  to  be  measured.  One 

advantage  of  periodic  excitation  methods  is  that  they  are  relatively  easily 

implemented.  Their  major  shortcoming  is  that  in  cases  where  several 
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Figure  1.  Three  basic  types  of  transonic  diffuser 
flows.  Mff  is  Mach  number  at  +  signs. 
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examine  the  mechanisms  responsible  for  the  unsteadiness  in  a  transonic 
diffuser  flow. 


waves  carried  by 
boundary  layers 

87-222-74 


Fig.  2  Waves  associated  with  self*execited 
transonic  diffuser  flow  oscillations. 
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2.  OBJECTIVES 


The  basic  concept  for  the  experimental  work  which  is  the  subject  of  this 
contract  is  shown  in  Fig.  3.  A  means  would  be  developed  by  which  an  abrupt 
pressure  rise  could  be  Imposed  and  subsequently  sustained  at  the  downstream 
end  of  a  transonic  diffuser  flow.  The  pulse  would  initiate  an  upstream-moving 
compression  wave  that,  upon  reaching  the  shock,  would  create  both  acoustic  and 
convective  reflections.  The  three  waves  would  then  be  separable  by  virtue  of 
their  individual  signatures  on,  for  instance,  wall-pressure  measurements,  as 
shown  in  Fig.  3- 

The  specific  tasks  for  this  contract  are  as  follows. 

Task  1  Develop  a  device  capable  of  generating  abrupt,  ramp-type  pressure 
pulses  in  a  supercritical  diffuser  flow.  The  device  will  be 
compatible  with  an  existing  MDRL-provided  diffuser  model. 


Fig.  3  Features  of  strong>shock  diffuser  flow  and 
conceptuai  description  of  pulse  propagation. 
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Task  2  Connect  the  pulse  generator  to  the  diffuser  model  and  install  in  test 
facilif".  Install  a  traversable,  fast-response,  total/static- 
pres'' ^e-probe  instrumentation  system  and  interface  to  the  MDRL 
iaooratory  computer.  Develop  the  specialized  control  and  software 
required.  Check  out  all  systems  and  conduct  preliminary  tests  to 
verify  correct  operation  of  the  computer  and  inlet  model  and  obtain 
preliminary  data. 

Task  3  Conduct  test  program:  acquire  and  reduce  final  data  for  entire  shock- 
strength  range.  Make  results  available  to  other  NWC-designated 
investigators  as  early  as  possible.  Document  results. 

The  development  of  the  technique  (Tasks  1  and  2)  has  been  documented  in 
the  open  literature  (Pubs.  1,  2)  and  only  a  brief  review  will  be  given  here. 
The  purpose  of  this  report  is  to  provide  an  overview  of  the  entire  program, 
assemble  the  results  of  the  various  task  phases  and  activities  into  a  single 
publication,  and  document  the  final  experimental  results  (Task  3). 
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3.  EXPERIMENT 


3.1  Diffuser  Model  and  Flowfield 

The  diffuser  model  used  in  this  study  is  a  nominally  two-dimensional, 
converging-diverging  channel  (Fig.  4).  The  channel  has  an  exit-to-throat  area 
ratio  of  1.52.  Side-  and  bottom-wall  boundary  layer  growth  is  limited  by 
three  sets  of  forward-facing  (ram)  suction  slots.  The  top-wall  boundary  layer 
is  tripped  6.92  throat -heights  upstream  of  the  throat.  The  model  side  walls 
are  fitted  with  schlieren-quality  glass  panels  that  allow  the  employment  of 
optical  diagnostic  techniques. 

Two  reference  locations  are  defined.  The  exit  station  at  x/h*  *  8.65, 
within  the  constant-area  section  of  the  channel  downstream  of  the  divergent 
section,  has  been  used  to  provide  the  downstream  boundary  condition  for 
computational  analyses  intended  to  describe  the  flow.  A  second  (pulse) 


Side  view 


Pulsed 


Top  view 


191  180  23  178  20.3  15.1  165 


Dimensions  in  mm. 

Vertical  dimensions  doubled. 

Slot  sizes  enlarged  for  clarity. 

Percentages  denote  area  decrease  at  slots. 

17-122-76 

Fig.  4  Diffuser  model. 
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reference  location  at  x/h*  »  15.2  is  the  station  at  which  the  initial 

conditions  of  the  injected  pulse  are  determined. 

Of  the  three  types  of  flow  shown  in  Fig.  1,  only  the  top  (weak-shock)  and 

bottom  (strong-shock)  cases  are  observed  in  this  diffuser.  Spark  schlieren 

photographs  of  the  experimental  flowfields  corresponding  to  these  two  cases 

are  displayed  in  Fig.  5.  The  weak-shock  case  (M  -  1.235)  has  thin  boundary 

0 

layers  and  attached,  well-behaved  flow.  The  strong-shock  case  (M  -  1.353) 

0 

displays  a  characteristic  massive  separation  with  the  top  and  bottom  wall 
boundary  layers  merging  after  approximately  five  throat-heights. 

In  addition  to  the  obvious  differences  between  the  mean  flowfield 
properties  of  these  two  flows,  there  are  also  significant  differences  in  their 
dynamic  characteristics. 

For  the  weak-shock  case,  the  power-spectral- density  (PSD)  distribution  of 

any  measured  quantity  shows  up  to  three  distinct  peaks  which  are  associated 

with  large-scale  coherent  oscillations  extending  throughout  the  entire 

1  4 

subsonic  flowfield.  The  frequencies  associated  with  these  modes  of 
oscillation  are  adequately  predicted  on  the  basis  of  simple  linear  acoustic 
theory. 

For  the  strong-shock  case,  the  PSD  contains  only  one  peak,  which  is  not 
describable  by  acoustic  theory.  The  conjecture  has  been  that  the  mode  of 
oscillations  for  this  case  is  comprised  of  several  elementary  waves,  which  may 


<7-222-77 


Fig.  5  Spark'Schlieren  photographs  of  attached 
and  separated  flow,  (a)  pg/p^s  1.225  (weak  shock), 
b)  Po/p,*  1.395  (strong  shock).  Grid  spacing  is 
24.5'nini  horizontally  and  12.7-nim  vertically. 
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include  upstream-  and  downs tream-traveling  acoustic  waves,  and  a  convective 

1 8 

boundary  layer  disturbance.  The  boundary  conditions  for  the  acoustic  waves 

are  an  open-ended  duct  at  the  downstream  end,  and  the  shock  at  the  upstream 

end.  The  acoustic  reflection  at  the  shock  is  determined  by  the  Rankine- 

Hugoniot  conditions,  as  they  apply  to  a  compression  wave  incident  upon  a 

normal  shock.  This  interaction  also  displaces  the  shock  upstream,  which 

introduces  a  corresponding  change  in  the  boundary-layer  properties.  These 

changes  are  convected  downstream  at  some  significantly  lower  speed  than  the 

core-flow  velocity.  Evidence  of  the  existence  of  this  convective  wave  has 

17 

been  extensively  documented. 

Additional  details  of  the  mean,  self-excited,  and  externally  excited 
flowfields  in  the  diffuser  are  presented  in  Refs.  14-17. 

3.2  Pulse  Generator 

The  pulse-generating  system  was  required  to  be  capable  of  generating  a 
ramp-type  pressure  disturbance,  i.e.,  a  transition  to  a  higher  pressure.  The 
magnitude  of  the  rise  must  be  controllable,  and  must  be  maintained  for  the 
length  of  time  required  for  the  disturbance  to  travel  upstream  and  for  any 
reflections  to  travel  downstream,  as  indicated  in  Fig.  3.  The  total  travel 
time  is  typically  on  the  order  of  5-10  ms.  The  rise  time  must  be  considerably 
shorter;  the  design  target  was  0.1  ms.  The  process  of  perturbing  the  flow 
must  be  repeatable  in  a  relatively  short  period  of  time,  since  several  pulses 
will  be  required  in  ensemble-averaging  the  transducer  signals  to  improve 
signal-to-noise  ratio  and  to  eliminate  pulse-to-pulse  variations. 

The  concept  adopted  for  the  pulse-generation  system  is  shown  in  Fig.  6. 

The  flowfield  perturbation  in  the  diffuser  is  generated  by  the  abrupt 
injection  of  a  mass  flow  into  the  flow  path.  The  mass  flow  is  injected 
through  16  identical  nozzles  located  in  single  spanwise  rows,  one  each  on  the 
top  and  bottom  of  the  Injector  head.  The  air  source  for  the  nozzles  is  a 
shock  tube,  operated  as  a  shock  tunnel.  This  apparatus  provides  an 
appropriately  fast  rise  time;  however,  the  duration  of  the  pulse  is  limited  to 
the  usable  test  time  of  the  shock  tunnel,  which  is  governed  by  the  shock  and 
expansion  wave  propagation  times  in  the  shock  tube.  To  lengthen  the  test 
time,  a  secondary,  blow-down  air  supply  systan  with  a  sliding-block-type  valve 
is  employed.  The  shock  activates  the  sliding-block  valves.  The  secondary 
systan  does  not  have  the  fast  rise-time  of  the  shock  tube,  but  is  capable  of 
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Secondary  air  supplied  through 
side  wall  sustains  blockage 
after  shock  process  completed 


■Primary  jets,  induced  by 
shock-compressed  air, 
block  flow  abruptly 


Diaphragm 
foil  material 
on  roll 

Pneumatic  driver  actuator 
moves  driver  to  allow  feed 


G^^hock  tube  driver 


Foil-feed 

mechanism  87-222-78 


Fig.  6  Pulse  generating  system  details. 


extremely  long  pulse  durations,  depending  on  the  secondary  air  source.  The 
secondary  air  is  injected  into  the  diffuser  flow  through  a  second  set  of 
spanwise  holes  on  the  top  and  bottom  of  the  injector  head  that  are  identical 
to,  and  just  upstream  of,  the  primary  injection  jets  from  the  shock  tube.  The 
two  systems  working  together  satisfy  the  requirements  of  fast-rise-time  and 
long-duration  pulses.  A  more  detailed  description  of  the  pulse  system  is 
provided  in  Pubs.  1  and  2.  Figure  7  is  a  photograph  of  the  pulse-injection 
system  attached  to  the  diffuser  model. 

The  pressure  ratio  in  the  diffuser  (ratio  of  stagnation  to  exit  static 
pressures,  Pq/Pq)  is  controlled  by  the  two  hinged  throttle  flaps,  one  for  each 
passage  around  the  pulse  injector  head.  The  flaps  are  also  used  to  compensate 
for  asymmetries  in  the  diffuser  flow.  Splitting  the  diffuser  flow  halves  the 
time  required  for  the  jet  fluid  to  traverse  the  diffuser  flow.  However, 
splitting  the  flow  can  also  create  a  problem  in  that  the  upstream  travel  time 
of  the  induced  pulse  in  the  top  and  bottom  channels  may  be  different  if  the 
two  channels  contain  different  flow  velocities.  Such  is  the  case  in  the 
strong-shock  flow  where  the  diffuser  flowfield  is  asymmetric  from  top  to 
bottom.  Care  was  taken  to  adjust  the  flaps  such  that  the  two  pulses  arrived 
at  the  pulse  reference  location  (x/h*  -  15.2)  simultaneously. 
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87-222-79 

Fig.  7  Diffuser  model  and  pulse  generating  system. 

Figure  8  shows  traces  from  transducers  mounted  on  the  top  and  bottom  walls 
of  the  diffuser  model  at  the  pulse  reference  station.  The  rise  time  of 
approximately  75  microseconds  meets  the  targeted  value.  The  close  agreanent 
between  top  and  bottom  wall  pressure  histories  indicates  nearly  simultaneous 
arrival  of  the  pulses,  i.e.,  a  single  planar  wave  front  propagating  upstream. 


P' 
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Fig.  S  Top  and  bottom-wall  static  pressure 
histories,  x/\  >  1S.2 
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The  existence  of  a  truly  planar  wave  front  was  further  confirmed  by 
introducing  a  third  static- pressure  transducer  in  midstream  at  the  reference 
location,  and  verifying  that  its  measurement  of  the  pressure  rise  corresponded 
in  time  and  amplitude  to  the  wall  measurements. 

The  operation  of  the  pulse-generation  system  is  entirely  computer- 
controlled.  A  schematic  of  the  various  actuation  controls  is  presented  in 
Fig.  9.  After  each  shot,  the  driver  section  of  the  shock  tube  is  retracted 
and  the  diaphragm  foil  material  is  advanced.  The  sliding-block  valves  for  the 
secondary  air  supply  are  reset.  The  driver  section  is  closed  and 
pressurization  begins.  The  driver  chamber  pressure  is  monitored  by  the 
computer,  and  when  the  desired  pressure  is  reached,  a  knife  blade  is  thrust 
through  the  foil.  The  entire  cycle  takes  less  than  10  seconds,  sc  a 
relatively  short  run  time  is  required  to  acquire  sufficient  data  for  ensemble 
averaging . 
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A  set  of  typical  operating  conditions  for  the  pulse-injection  system  is 
provided  in  Table  1 .  Driver  pressures  from  70  kPa  to  350  kPa  were  generally 
used  for  the  experiments.  These  pressures  generated  pulse  strengths  which 
varied  over  a  3:1  range,  as  measured  at  the  pulse  reference  station. 


Table  1.  Typical  Shock  Tube  Operating  Conditions 


Pressure  in  driver  before  ruing 

375  kPa 

Pressure  in  tube  before  firing 

100  kPa 

Temperatures  in  tube  and  driver 

300  K 

Shock  speed 

460  m/s 

Pressure  of  doubly-shocked  gas 

333  kPa 

Temperature  of  doubly-shocked  gas 

428  K 

Jet  Mach  number 

i.43 

Jet  speed 

500  m/s 

Duration  of  steady  jet  flow 

5.2  ms 

Time  for  jets  to  traverse  a 

channel  branch  (set  to  2.0-cm  height) 

40  (IS 

87.^2-71 


3«3  Instrumentation  and  Data  Acquisition 

The  diffuser  operating  condition  was  routinely  monitored  and  recorded, 
including  measurements  of  the  diffuser  supply  plenum  static  pressure  and 
temperature  (diffuser  total  pressure  and  total  temperature),  the  diffuser 
pressure  ratio,  and  atmospheric  pressure.  In  addition,  the  burst  pressure  of 
the  shock  tube  driver  was  recorded  for  each  shot  in  order  to  determine  the 
amount  of  shot-to-shot  variation.  With  the  computer-triggered  diaphragm  burst 
procedure,  this  variation  was  on  the  order  of  3.5  kPa  (rms). 

The  top  wall  of  the  diffuser  model  is  fitted  with  seven  flush-mounted, 
fast-response  pressure  transducers,  distributed  from  the  vicinity  of  the  shock 
to  just  upstream  of  the  injector  head.  The  pulse  reference  signals  are  used 
to  determine  the  pulse  strength,  which  is  used  in  the  data  reduction  to 
normalize  the  rest  of  the  wall  pressures.  In  addition,  a  fast-response 
transducer  was  located  in  the  driven  side  of  the  shock  tube  near  the  injector 
head  to  trigger  the  high-speed  data  acquisition  system. 

A  traversable  twin  total/static  "core"  probe  (similar  to  that  described  in 
Ref.  16)  fitted  with  fast-response  transducers  was  used  to  map  the  unsteady 
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pressure  field  in  the  core  flow.  The  probe  could  be  moved  vertically  along 

the  full  height  of  the  channel,  and  in  the  streamwise  direction  from  the  shock 

to  just  upstream  of  the  exit  station. 

Optical  data  were  also  acquired.  High-speed  schlieren  movies  (10,000 

frames/s)  were  made  to  document  the  shock  response  to  the  pulse  impingement. 

The  movies  became  a  valuable  diagnostic  tool,  especially  for  the  weak-shock 

case  in  which  it  was  found  that  the  same  pulse  strength  used  for  the  strong- 

shock  case  would  drive  the  weak  shock  completely  through  the  throat,  i.e.,  the 

pressure  rise  imposed  at  the  downstream  end  of  the  diffuser  generated  a 

pressure  ratio  corresponding  to  fully  subsonic  flow  in  the  diffuser 

(subcritical  operation).  Hence,  a  weaker  pulse  was  used  for  the  weak-shock 

flow  to  maintain  fully  supercritical  operation.  In  addition,  a  real-time 

24  25 

shock-position-indicating  system  ’  was  used  to  measure  the  time  history  of 
the  shock  location,  corresponding  to  the  time  histories  of  the  wall  and  core 
pressures . 

All  unsteady  data  were  digitized  at  50  kHz/channel.  This  provided  20- 
microsecond  resolution  which  was  deemed  sufficient.  Trials  with  digitization 
rates  up  to  100  kHz/channel  were  made;  no  significant  difference  appeared  in 
the  data.  Ensemble  averages  were  made  over  25  shots  for  all  data.  Trials 
were  made  for  ensemble  averaging  of  10  to  50  shots.  It  was  found  that  for 
more  than  25  shots,  no  significant  changes  or  smoothing  appeared  in  the  data. 

3.4  Test  Conditions 

The  shock  Mach-number  range  examined  was  1.18  <  M  <  1.38,  with  the 

0 

transition  from  weak-shock  to  strong-shock  flow  occurring  at  about  = 

1.273.  This  Mach-number  range  corresponds  to  a  pressure-ratio  range  of  1.21  < 

p^/Pg  <  1.41.  Wall-pressure  measurements  were  taken  at  10  pressure  ratios 

covering  tnis  entire  range.  Two  shock  Mach  numbers  (M  -  1.251  and  M  - 

0  0 

1.378)  were  selected  for  detailed  examination.  For  these  two  cases,  a  study 
was  made  of  the  effect  of  pulse  strength  on  the  pulse  propagation 
characteristics,  and  core-probe  surveys  were  taken. 

For  the  strong-shock  data,  the  driver  pressure  was  set  at  275  kPa.  This 
created  a  pulse  whose  strength  was  6.3%  of  the  local  time-mean  static  pressure 
as  measured  at  the  injector  head  reference  station.  For  the  weak-shock  core¬ 
probe  data,  the  driver  pressure  was  137  kPa,  creating  a  pulse  strength  of 
3.0$. 
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The  initial  runs  revealed  that  the  primary  pulse  provided  by  the  shock 
tube  system  had  a  sufficiently  long  duration  for  the  upstream  and  downstream 
waves  to  complete  their  travel.  The  secondary  air-injection  system  thus  was 
not  required  and  its  use  was  discontinued. 
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4.  RESULTS 


4.1  Wall-Pressure  Measurements 

A  sample  set  of  time  histories  of  the  top- wall  transducers  is  shown  in 
Fig.  10  for  a  weak-  (M  -  1.251)  and  strong-shock  flow  (M  -  1.378).  Passage 

<J  0 

of  the  upstream-traveling  acoustic  wave  can  clearly  be  seen  as  the  sharp  rise 
in  the  signal.  The  rise  occurs  progressively  later  at  each  successive 
upstream  measurement  station.  The  pulse  is  initially  quite  steep,  but 
broadens  as  it  approaches  the  shock. 

In  all  of  the  weak-shock  cases,  only  the  upstream- traveling  acoustic  wave 

appears;  there  is  no  observed  reflection  of  any  sort.  This  is  consistent  with 

linear  acoustic  theory,  which  predicts  negligible  acoustic  reflection  for  weak 
26 

shocks.  In  the  strong-shock  case,  however,  the  reflection  is  strong  and  can 
be  seen  clearly  as  a  secondary  rise  in  the  pressure  signals  (arrows)  ,  which 
occurs  at  later  times  at  the  more  downstream  stations. 

From  transducer-signal  time  histories  such  as  those  presented  in  Fig.  10, 
x-t  diagsrams  can  be  generated  from  which  the  wave  propagation  speeds  can  be 
determined.  Figure  11  contains  an  x-t  diagram  for  both  the  pulse  propagation 
(Fig.  11b)  and  the  response  of  the  shock  to  the  impingement  of  the  acoustic 
pulse  (Fig.  11a). 

Figure  1 1b  contains  the  pulse-propagation  information  for  a  weak  shock  (M^ 

=  1.251),  a  strong  shock  (M^  -  1.378),  and  a  case  where  the  shock  is  just 

strong  enough  to  cause  a  shock-induced  separation  (M  -  1.292).  It  was  found 

0 

that  the  speed  of  all  waves  was  nearly  constant  over  the  length  of  their 
travel;  the  x-t  diagrams  are  essentially  straight  lines.  The  reflected  waves 
in  the  strong-shock  cases  are  clearly  convective,  as  indicated  by  the 
relatively  slow  speed  of  95  m/s  (an  acoustic  wave  would  travel  4-5  times  as 
fast).  No  evidence  of  a  reflected  acoustic  wave  was  observed  for  any  of  the 
flows . 

If  the  strong-shock  case  acoustic-  and  convective- wave  trajectories  in 

Fig.  11  are  extrapolated  upstream,  then  the  locations  where  they  cross  are 

nearly  midway  between  their  respective  undisturbed  and  final  shock  locations. 

That  the  convective  wave  for  the  M  ■  1.292  case  is  delayed  somewhat  from  that 

0 

of  the  M  -  1.378  case  results  from  the  fact  that  the  undisturbed  shock 
0 

location  for  the  weaker  of  the  two  shocks  is  farther  upstream  and  the  acoustic 
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(a)  (b) 


Fig.  11  x*t  diagram  of  weak  and  strong-shock  flow  response  to 
pulse:  (a)  instantaneous  shock  location  (horizontal  scale 
expanded);  (h)  pulse  propagation  measured  along  top  wall. 

wave  must  travel  that  much  farther  before  encountering  the  shock.  This  was 
found  to  be  the  case  for  all  strong-shock  flows  examined.  There  is  no 
measurable  time  delay  between  the  impingement  of  the  acoustic  wave  and  the 
generation  of  the  convective  wave. 

Table  2  contains  a  summary  of  the  wave  propagation  speeds  for  all  the 
flows  examined.  In  spite  of  the  complex,  decelerating,  scmetimes  separated 
flowfields  in  each  of  these  cases,  the  propagation  speeds  seem  remarkably 
independent  of  pressure-ratio  effects.  The  upstream  speed  (w_)  is 
consistently  around  200  m/s,  and  the  downstream  speed  (w^)  is  approximately 
half  that.  The  observed  w_  values  are  unexpectedly  high,  especially  near  the 
shock  where  the  core-flow  velocity  is  the  highest.  This  phenomenon  is 
examined  in  more  detail  in  Section  5.2. 

The  shock  response  to  the  impinging  acoustic  pulse  (Fig.  11a)  shows 
different  characteristics  for  the  weak-  and  strong-shock  cases.  For  the  weak 
shock,  the  shock  responds  almost  instantaneously  as  the  acoustic  pulse 
impacts,  but  the  time  required  for  the  shock  to  assume  its  new  equilibrium 
location  corresponding  to  the  lower  pressure  ratio  created  by  the  pulse  is 
relatively  long.  On  the  other  hand,  for  the  strong  shock,  there  is  a  delay  of 
1  -  1.5  ms  before  the  shock  responds  to  the  acoustic  wave,  by  which  time  the 
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Table  2.  Wave  Propagation  Speeds 


Ma 

Po^Pe 

(m/s) 

(m/s) 

Strong  Shock: 

1.378 

1.410 

199 

90 

1.339 

1.383 

200 

90 

1.320 

1.359 

192 

100 

1.309 

1.335 

200 

96 

1.292 

1.309 

199 

102 

1.276 

1.280 

202 

117 

Weak  Shock: 

1.270 

1.258 

205 

— 

1.251 

1.232 

205 

— 

1.177 

1.207 

202 

— 

S7.222.72 


reflected  convective  wave  has  traveled  nearly  five  throat-heights  downstream. 
However,  once  the  shock  does  respond,  it  assumes  its  new  equilibrium  location 
quite  quickly. 

A  study  was  also  made  of  the  effect  of  pulse  strength  variation  for  the 
cases  of  Fig.  10.  The  initial  pulse  strength,  defined  as  the  ratio  of  the 
pressure-rise  magnitude  at  the  pulse  reference  station  to  the  mean  pressure  at 
the  same  location,  was  varied  over  a  more  than  3:1  range.  The  results  are 
shown  in  Fig.  12.  Normalizing  the  pulses  measured  along  the  length  of  the 
model  (p')  with  the  pulse  reference  station  strength  (p'  produces  a  good 

X  I  b .  t 

correlation  of  the  data  for  both  cases. 

For  the  weak-shock  case  (Fig.  12a),  the  pulse  is  somewhat  attenuated  as  it 
travels  upstream,  but  rises  sharply  as  the  shock  is  approached.  Similar 
trends  were  also  found  when  a  periodic  excitation  was  applied  to  the 
downstream  end  of  the  same  diffuser  model  (Figs.  9-11  of  Ref.  1 6) .  This 
similarity  between  the  pulsed  and  periodic  excitation  is  not  surprising  since 
for  weak-shock  flows  only  the  upstream-traveling  acoustic  wave  is  present, 
whether  the  perturbation  is  pulsed  or  periodic. 

For  the  strong-shock  case,  the  upstream  wave  is  slightly  more  attenuated 
than  for  the  weak  shock.  The  reflected  convective-wave  amplitude  is  large 
Initially,  but  strongly  damped,  diminishing  to  zero  at  approximately  x/h*  - 
11.  This  location  is  near  several  significant  mean-flow  features.  The  jore 
flow  terminates  at  x/h*  -  8.4,  and  the  separation  bubble  terminates  at  x/h*  - 
6.1.  Neither  of  these  two  features  exists  in  the  weak-shock  flow,  which 
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Fig,  12  Normalized  pulse  strength  distributions: 

(a)  1.252;  (b)  1,378.  Open  symbols: 

upstream-moving  acoustic  waves;  full  symbols: 
downstream-moving  convective  waves:  balf-full 
symbols:  contributions  from  both  waves  present. 

suggests  that  a  thlck-boundary-layer/core-f low  interface  and  a  shock- induced 
separation  may  both  be  related  to  the  existence  of  the  convective  wave.  In 
addition,  from  the  shock  to  approximately  x/h*  -  n,  the  velocity  profiles 
display  inflection  points  in  the  top-wall  boundary  layer.  This  particular 
feature,  for  a  free  shear  layer  between  streams  of  infinite  extent,  is  known 
to  ensure  the  growth  of  disturbances  over  a  finite  frequency  range.  In  the 
present  situation  (wall  proximity,  decelerating  flow,  and  increasing 
transverse  spatial  scale)  the  disturbance  does  not  grow,  but  the  damping  is 
moderate  enough  to  allow  observation  of  the  disturbance. 

4.2  Core-Flow  Pressure  Measurements 

Both  total-  and  static-pressure  measurements  were  made  throughout  the  flow 
from  just  downstream  of  the  shock  to  Just  upstream  of  the  exit  station  in  both 
the  weak-  and  strong-shock  flowfields.  The  total-pressure  measurements  are 
difficult  to  interpret  alone  since  they  are  composed  of  fluctuations  in  both 
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static  pressure  and  velocity.  Since  detailed  analysis  is  beyond  the  scope  of 
this  contract,  only  the  static-pressure  data  will  be  discussed. 

For  the  weak-shock  flow,  the  core-probe  static- pressure  signals  were 

similar  to  the  wall-pressure  signals.  The  core-probe  signals  showed  only  a 

single  abrupt  rise  corresponding  to  the  upstream- traveling  acoustic  wave. 

Data  from  the  survey  have  been  assembled  to  display  the  propagation  of  the 

wave  front.  The  wave  front  was  said  to  have  arrived  when  the  signal 

significantly  exceeded  the  pre-pulse  noise  level  (about  1  kPa) .  Figure  13 

shows  the  wave  front  location  at  every  0.1  ms.  The  wave  front  remains 

essentially  vertical  throughout  its  travel,  indicating  that  the  pulse  is 

1 5 

plana”,  despite  the  nonuniforra  velocity  profile.  The  wave  speed  as 

determined  from  the  core-probe  data  is  the  same  as  that  determined  from  the 
wall  data. 

The  strong-shock  core-probe  data  do  not  provide  as  simple  a  picture  as  the 
weak-shock  data.  Figure  14  shows  the  core-probe  static-pressure  signals  at 
x/h*  =>  4.0  at  selected  locations  spanning  the  height  of  the  channel.  Near  the 
top  wall  the  core-probe  data  closely  resemble  the  wall  data;  however, 
additional  ripples  appear  in  the  signal,  becoming  more  sharply  defined  as  the 
bottom  wall  is  approached.  The  ripples  display  a  well-defined  pattern  in 
streamwise  variation  as  well.  Figure  15  shows  core-probe  pressure  traces  for 
a  streamwise  survey  at  y/h*  »  0.3.  As  the  shock  is  approached,  the  number  of 
peaks  which  appear  increases  from  one  at  the  downstream  end  of  the  surveyed 


y/h^  87-22245 

Fig.  13  Weak'Shock  pulse  wavefront  propagation.  Contours  at  0.1  ms  intervals. 
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Fig.  14  Core  static  pressure  measurements,  verticai  scan  at 
x/h,  =  4.0. 


region  to  three  near  the  shock.  The  spacing  between  the  peaks  changes  as 
well,  indicating  that  the  ripples  are  dispersive. 

An  attempt  was  made  to  determine  the  origin  of  the  ripples.  Probe 
interference  or  mechanical  vibration  effects  were  investigated  and  ruled  out 
as  the  cause.  An  x-t  diagram  (Fig.  l6a)  was  constructed  containing  the 
initial  pulse  arrival  (determined  as  for  the  weak-shock  case  described  above) 
and  the  first  and  second  peaks,  along  with  the  convective  wave.  A  least- 
squares  fit  was  made  to  the  data.  The  speed  associated  with  the  pulse  arrival 
(203  m/s)  is  essentially  the  same  as  that  for  the  pulse  as  measured  on  the  top 
wall,  while  the  speeds  associated  with  the  first  and  second  peaks  are  135  m/s 
and  93  m/s ,  respectively.  Throughout  the  range  where  a  third  peak  could  be 
clearly  discerned,  its  speed  was  measured  to  be  79  m/s. 
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Fig.  15  Core  static  pressure  disturbance 
propagation,  horizontal  scan  at  y/h.  =  0.3. 

A  striking  feature  of  the  peak  trajectories  is  that  they  all  seem  to 
converge  at  a  single  point,  approximately  x/h*  -  9.6l  .  Figure  l6b  shows 
various  features  of  the  flowfield  and  model  geometry.  A  suction  slot  is 
located  at  x/h*  -  9.80,  quite  close  to  the  intersection  of  the  trajectories. 
The  possibility  that  the  suction  slot  and  its  associated  plenum  were  acting  as 
a  Helmholtz  resonator  was  tested  and  ruled  out,  as  was  the  possibility  that 
the  ripples  originated  downstream  at  the  injector  head  and  the  convergence  of 
the  trajectories  at  the  slot  is  Just  a  coincidence.  The  final  conclusion  is 
that  the  abrupt  12..7%  area  change  at  the  bottom  slot  triggers  the  ripples; 
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<«t» 


Fig.  16  Tracing  origin  of  ripples:  (a)  x>t  diagram 
of  propagation;  (b)  modei  and  flow  features. 


their  prominence  near  the  bottom  wall  is  probably  related  to  the  fact  that  the 
suction  slots  exist  only  on  the  bottom  and  side  walls,  not  on  the  top  wall. 

The  upstream-moving  compression  wave,  as  it  exists  for  x/h*  >  9.80,  is  the 
intended  perturbation.  The  passage  of  the  pulse  through  the  abrupt  area- 
change  at  the  slot  constitutes  an  additional  perturbation  modifying  the  pre¬ 
existing,  intended  perturbation.  This  two-stage  triggering  process  is 
complex,  as  is  the  resulting  wave-front  structure.  The  evolution  of  the 
composite  perturbation  as  evidenced  by  the  ripples  on  the  pressure  signals  is 
reminiscent  of  nonlinear  gravity  waves  on  liquid  surfaces. 

The  mechanism  for  triggering  the  ripples  is  model-specif ic,  but  their 
positive  growth  rate  suggests  that  they  might  occur  in  other  situations  as  a 
result  of  other  triggering  mechanisms.  The  nature  of  the  ripples  is  not 
understood,  but  further  Investigation  of  this  phenomenon  is  beyond  the  scope 
of  this  contract. 

The  ripples  obscure  the  pressure  signature  of  the  convective  wave  over  a 
significant  portion  of  the  channel.  Figure  I6a  is  derived  from  data  taken  at 
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y/h*  =  0.5.  In  these  data  the  convective  wave  cannot  be  observed  until  it  is 
nearly  four  throat-heights  downstream  of  the  shock.  Farther  upstream,  the 
convective  wave  is  overlapped  by  the  oscillations.  As  seen  in  Fig.  14,  the 
anomalous  peaks  decrease  in  amplitude  toward  the  top  wall,  but  the  actual 
arrival  of  the  convective  wave  cannot  be  determined  with  confidence  until  very 
near  the  top  wall.  However,  at  the  measurement  stations  nearest  the  top  wall, 
the  wave  speeds  measured  by  the  core  probe  correspond  closely  to  those 
obtained  frcm  the  wall-pressure  measurements. 
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5.  RELATED  ANALYSIS 


Inspired  by  the  experimental  results  obtained  through  this  contract  work, 
we  performed  several  related  studies  under  funding  by  the  McDonnell  Douglas 
Independent  Research  and  Development  program.  These  studies  extended  the  data 
analysis  beyond  the  scope  of  the  contract,  and  addressed  some  appeirently 
anomalous  behavior  in  the  results.  These  studies  are  briefly  described  below. 

5.1  Resultant  Wave  Reconstruction 

Using  the  amplitude  and  phase  information  from  the  experimental  data 
obtained  under  this  contract,  an  attempt  was  made  to  reconstruct,  from  the 
elementary  waves,  the  oscillations  found  in  the  strong-shock  case  where  the 
flow  was  periodically  excited.  (The  weak-shock  case  is  trivial,  since  only 
one  wave  is  present.)  Experimental  data  from  Ref.  16,  where  the  flow  was 
excited  at  100,  200,  and  300  Hz,  are  used  for  comparison.  The  success  of  such 
a  reconstruction  would  prove  that  the  elementary  wave  properties  determined 
experimentally  from  the  pulse  tests  form  a  valid  basis  for  modeling  a  variety 
of  unsteady  flows. 

The  reconstruction  was  accomplished  by  forming  two,  counter propagating, 
simple,  single-frequency  waves  having  the  streamwise  amplitude  distributions 
and  speeds  measured  for  the  elementary  waves  from  this  pulsed  experiment.  The 
individual  wave  numbers  were  derived  from  the  wave  speed  and  frequency.  These 
two  waves  were  then  simply  added  together.  Special  considerations  involved  in 
determining  the  relative  amplitudes  and  phases  of  the  component  waves  are 
discussed  below. 

Separation  of  amplitude  contributions  from  the  upstream-  and  downstream¬ 
traveling  waves  on  the  basis  of  the  wall-pressure  measurements  alone  is  not 
possible  near  the  shock.  The  time  delay  between  the  two  waves  as  they  pass 
over  the  transducer  is  negligible  and  the  signal  shows  only  a  single  rise 
(Fig.  10b,  x/h*  -  3-13).  The  uncertainty  of  the  relative  contributions  from 
the  two  waves  at  this  location  is  emphasized  by  the  half-filled  symbols  in 
Fig.  12b.  The  convective- to- acoustic  amplitude  ratio  at  the  shock,  based  on 
the  data  trends  farther  downstream,  was  estimated  to  be  3j2.  (Further 
analysis  of  the  core-probe  data  is  required  to  refine  this  estimate.) 
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The  phase  relation  between  the  two  waves  was  initially  determined  by 
considering  the  x-t  diagram  in  Fig.  11b.  The  negligible  time  delay  between 
the  impingement  of  the  acoustic  wave  and  generation  of  the  convective  wave 
indicates  that  there  is  no  phase  shift  between  the  upstream-  and  downstream¬ 
traveling  waves  at  the  shock.  This  observation,  coupled  with  the  estimated 
amplitude  relation  discussed  above,  yields  a  reflection  coefficient  of  1.5; 
and  since  the  pulse  reflection  appears  with  no  time  delay  and  has  the  same 
sense  as  the  incident  wave  (pressure  rise),  the  reflection  coefficient  for  a 
harmonic  wave  based  on  the  pulse  response  would  be  a  real,  positive  number. 

It  is  not  certain  to  what  degree  the  phase  relations  implied  in  the  pulse 

experiment  can  be  extended  to  the  case  of  harmonic  waves;  the  spectral 

contents  of  the  two  waves  are  completely  different-  Analysis  of  normal  shock 
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response  to  harmonic  acoustic  waves  predicts  a  frequency-dependent  phase 
shift  at  the  shock.  The  considerable  delay  between  the  arrival  of  the 
acoustic  wave  and  the  eventual  recoil  of  the  shock  (Fig.  11)  also  implies  that 
some  additional  phase  shift  may  occur  in  the  harmonic  case.  To  optimize  the 
agreement  between  the  reconstruction  and  the  experimental  data,  an  additional 
delay  of  0.5  ms  between  the  upstream-  and  downstream-traveling  waves  at  the 
shock  was  incorporated  in  the  reconstruction  for  all  cases. 

The  results  of  the  wave  reconstructions  ah's  shown  in  Fig.  17.  The 
experimental  phase  data  were  referenced  to  the  angular  position  of  the 
mechanical  excitation  device;  hence,  the  relative  phase  relation  between  the 
data  and  the  reconstruction  is  unknown  since  no  equivalent  reference  exists 
for  the  pulse  experiment.  The  reconstructed  phase  data  were  arbitrarily 
shifted  to  agree  with  the  experimental  data  near  the  shock,  with  the  same 
shift  applied  to  the  rest  of  the  downstream  data.  In  addition,  it  has  been 
the  practice  with  this  model  to  normalize  amplitude  data  with  the  amplitude  at 
the  exit  station.  Figure  17  follows  this  practice;  the  experimental 
amplitudes  and  the  reconstruction  have  been  normalized  with  their  respective 
values  at  x/h*  -  8.65. 

The  agreement  of  the  phases  between  the  experimental  data  and  the 
reconstruction  is  excellent  for  the  100-Hz  and  200-Hz  cases.  The  agreement 
for  the  300-Hz  case  is  less  striking,  although  the  major  trends  in  the 
reconstruction  still  follow  those  of  the  experimental  data.  The  major 
discrepancy  is  the  failure  of  the  reconstruction  to  capture  the  sharp  phase- 
shift  at  the  node  at  x/h*  -  3.8.  However,  a  sharp  phase-shift  can  only  occur 
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if  there  is  exact  cancellation  between  the  two  counter-propagating  waves, 
which  requires  that  the  two  waves  have  identical  amplitudes  at  that  point. 
Given  the  uncertainty  in  estimating  the  relative  wave  amplitudes  near  the 
shock,  it  is  not  surprising  that  this  particular  feature  is  not  well 
reproduced . 

A  heretofore  puzzling  feature  of  the  experimentally  determined  periodic 
phase  distributions  is  the  evidence  of  an  upstream  propagating  wave 
(decreasing  phase  with  increasing  distance)  past  the  end  of  the  core  flow 
(x/h*  >  8)  present  in  all  three  cases.  It  was  expected  that  the  externally 
imposed  periodic  excitation  would  generate  a  simple  standing- wave-1 ike 

1 8 

amplitude  envelope  (nodes  and  antinodes)  superimposed  over  a  traveling  wave. 
However,  given  the  results  shown  in  Fig.  12b,  the  reason  for  the  observed 
behavior  becomes  evident.  Downstream  of  the  core  flow,  the  downstream¬ 
traveling  convective  wave  has  completely  died  out.  The  ups tream- travel ing 
acoustic  wave  is  the  only  one  present  and  hence,  the  only  influence  observed. 

The  amplitude  distributions  for  the  experimental  data  and  the 
reconstructions  also  show  qualitative  agreement;  the  trends  of  maxima  and 
minima  are  faithfully  represented;  however,  there  is  considerable  difference 
in  the  relative  magnitudes.  For  the  100-Hz  case,  the  reconstructed  amplitudes 
are  large  everywhere.  Near  the  shock,  the  reconstructed  peak  is  nearly  an 
order  of  magnitude  larger  than  the  data.  The  300-Hz  case  shows  only  about  a 
factor  of  two  difference  between  the  data  and  the  reconstruction.  At 
x/h*  »  6,  the  reconstructed  peak  is  larger,  but  the  apparent  approach  to  a 
peak  just  past  x/h*  -  11  shows  the  data  to  be  larger.  The  200-Hz  case  shows 
excellent  agreement.  (It  might  be  noted  that  200  Hz  is  near  the  self-excited 
oscillation  frequency  for  this  diffuser,  217  Hz.) 

Finally,  with  regard  to  the  self-excited  oscillations,  Ref.  14  found  that 

for  the  strong-shock  flow  the  self-excited  frequencies  could  not  be  described 

by  simple  acoustic  theory.  The  frequencies  did  not  scale  with  the  length  of 

the  subsonic  flowfield,  but  rather  with  the  core-flow  length.  Using  the  wave 

speeds  measured  in  this  pulsed  experiment,  predictions  were  made  of 

frequencies  by  calculating  the  round-trip  travel  time  based  on  the  measured 
1  4 

core-flow  length  .  This  procedure  implies  that  the  convective  wave  generates 
an  upstream-traveling  acoustic  wave  at  the  downstream  end  of  the  core  flow  by 
some  '’as  yet  unknown)  mechanism.  The  acoustic  wave,  in  turn,  generates 
another  convective  wave  when  it  reaches  the  shock.  The  results  are  shown  in 
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Fig.  1 8.  The  agreement  between  the  self-excited  data  and  the  prediction  is 
quite  good,  strongly  suggesting  that  the  mechanism  responsible  for  the  self- 
excited  oscillations  consists  of  a  cycle  involving  the  convective  and  acoustic 
waves . 

This  analysis  of  the  pulsed  data  shows  that  the  measured  oscillations 
occurring  in  both  the  self-excited  and  externally  excited  strong-shock 
diffuser  flow  result  from  the  interaction  of  elementary  waves,  a  fact  which 
until  now  has  only  been  conjecture.  It  also  indicates  that  these  elementary 
waves  are  the  ups tr earn- travel ing  acoustic  wave  and  a  downs tr earn- travel ing 
convective  wave  observed  in  the  pulsed-flow  measurements. 

5.2  Pulse  Propagation  in  Nonunifom  Channel  Flows 

An  anomalous  finding  described  in  section  4.1  is  that  the  speed  of  the 
upstream-propagating  acoustic  wave  is  considerably  faster,  especially  near  the 
shock,  than  would  be  predicted  on  the  basis  of  core-flow  properties.  Since 
knowledge  of  the  wave  speed  is  necessary  to  calculate  the  self-excited 
frequencies  of  a  diffuser,  a  simple  analytical  method  for  predicting  the  wave 
speed  was  developed  for  nonuniform  channel  flows. 


87-222-90 

Fig.  18  Comparison  of  self-excited  frequency 
prediction  based  on  elementary  wave  properties 
obtained  from  experimental  data. 
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Viscous  effects  enter  into  the  analysis  only  in  that  they  are  responsible 
for  the  nonuniformity  of  the  velocity  and  density  profiles;  they  are  otherwise 
neglected.  The  pressure  is  assumed  constant  normal  to  the  flow,  and  the  shape 
of  the  cross  section  and  the  flow  property  profiles  may  be  specified 
arbitrari ly . 

Weak  compression  waves  are  treated  by  requiring  the  conservation  of  mass, 
momentum,  and  energy,  all  of  which  are  expressed  as  appropriate  integrals  over 
the  entry  and  exit  surfaces  of  a  control  volume  fixed  to  the  wave.  The 
integrals  are  expressed  in  terms  of  shape  factors  and  mean-flow  properties. 

The  major  assumption  made  is  that  the  shapes  of  the  density  and  velocity 
profiles  (and  hence  the  shape  factors)  are  the  same  on  both  faces  of  the 
control  volume;  only  the  magnitude  of  these  quantities  is  changed  by  the 
passage  of  the  wave.  Since  the  shape  factors  are  known  from  the  specification 
of  the  initial  profiles,  the  number  of  unknown  properties  behind  the  wave  is 
reduced  to  three,  and  can  thus  be  computed  from  the  three  available  equations. 

Through  the  above  formalism,  the  wave  speed  can  be  determined  by  use  of 
shape  facto''s  computed  from  the  velocity  profiles  measured  relative  to  the 
wave.  Since  the  wave  speed  is  not  known  initially,  the  wave  speed  and  the 
shape  factors  appropriate  for  the  coordinate  system  moving  with  the  wave  are 
determined  simultaineously . 

Results  of  applying  the  above  technique  to  several  numerical  test  cases 
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show  good  agreement  with  numerical  predictions.  Figiore  19  shows  the 
dependence  of  wave  speed  on  boundary-layer  thickness  for  a  simple  profile 
constructed  of  a  constant-speed  core  flow  and  a  parabolic  boundary-layer. 

This  technique  was  also  applied  to  the  weak-  and  strong-shock  flowfields 
examined  in  this  study.  Figure  20  shows  a  comparison  of  the  experimental 
data,  the  wave  speed  predicted  through  application  of  1-D  gas  dynamics  using 
the  core  flow  conditions,  and  the  prediction  based  on  the  present  theory. 
Although  a  minor  discrepancy  still  remains  near  the  shock,  the  present  theory 
predicts  the  wave  speed  remarkably  well.  The  comparison  Indicates  that 
although  the  theory  was  derived  for  constant-area  channels,  it  is  adequate  for 
cases  in  which  the  cross-sectional  area  varies  slowly,  provided  that  the 
calculation  is  based  on  local  flow  profiles.  Details  of  this  theory  are 
presented  in  Ref.  28. 
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Fig.  19  Dependence  of  normalized  wave  speed  on 
boundary  layer  thickness.  Core  flow  Mach  number 
is  0.7.  (From  Ref.  28). 


NWC  TP  6846 


Q  ■  Data 

•  ”  “  “  Expectation  based  on  core-flow  speed 
^  ■  Present  theory 


87-222-92 


Fig.  20  Trajectory  of  opstream-moving  wave. 
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6.  SUMMARY 

The  significant  accomplishments  and  results  of  this  program  are  the 
following. 

»  A  fast- rise- time  pulse  generator  was  designed  and  constructed,  and 
operated  according  to  expectations. 

*  When  used  with  a  transonic  diffuser,  the  pulse  generation  technique  was 
capable  of  separating  the  elementary  waves  for  individual  investigation. 

*  Of  the  three  elementary  waves  conjectured  to  exist  in  transonic  diffuser 
flows,  only  two  (upstream- traveling  acoustic  and  downstream- traveling 
convective  waves)  were  found  to  exist  in  the  strong-shock  flow,  while  only 
the  upstream- traveling  acoustic  wave  was  found  in  the  weak-shock  flow. 

»  The  amplitude  distributions  and  phase  relations  of  the  elementary  waves 
were  determined. 

*  The  pulse  reflection  coefficient  at  the  shock  for  the  weak-shock  case  was 
found  to  be  zero,  while  for  the  strong-shock  case  it  was  estimated  to  be 

1.5. 

*  The  convective  wave  was  strongly  damped  and  existed  only  in  those  regions 
where  the  velocity  profile  possessed  an  inflection  point. 

*  From  the  elementary  wave  information  determined  in  this  program,  features 
of  the  resultant-  wave  patterns  for  periodically  excited  strong-shock  flow 
could  be  reconstructed  with  reasonable  agreement,  and  the  frequencies  for 

I  the  self-excited  strong-shock  flow  could  be  predicted. 
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Core-flow  measurements  for  the  weak-shock  case  showed  a  planar  (one- 
dimensional)  acoustic  wavefront  structure  at  all  locations  in  the  flow, 
with  wave  speeds  corresponding  closely  to  those  obtained  fran  wall- 
pressure  measurements. 

Core-flow  measurements  for  the  strong-shock  case  revealed  a  family  of 
rapidly  growing  ripples  superimposed  on  the  original  pulse. 
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7.  PRESENTATIONS  AND  PUBLICATIONS 

(based  completely  or  in  part  on  work  done  under  this  contract) 

1.  Sajben,  M.  and  Kroutil,  J.  C.,  "Concept  and  Technique  for  Probing  Channel 
Flows  with  Abrupt  Perturbations”,  AIAA  Paper  No.  86-0311.  January  1986. 

2.  A  revised  version  of  the  above  paper  with  the  same  title  has  been  accepted 
by  the  Journal  of  Propulsion  and  Power. 

3.  Sajben,  M.,  Bogar ,  T.  J. ,  and  Atkins,  H.  L. ,  "Propagation  of  Weak  Shocks 
in  Nonuniform  Channel  flows",  US/French  Canbustion  Instability  Workshop, 
Naval  Weapons  Center,  China  Lake,  CA,  April  1986. 

4.  Bogar,  T.  J.  and  Sajben,  M.,  "Response  of  Transonic  Diffuser  Flows  to 
Abrupt  Increases  of  Back  Pressure;  Wall  Pressure  Measurements,  23rd 
JANNAF  Combustion  Meeting,  October  1986,  CPIA  Pub.  457. 

5.  Because  of  the  limited  distribution  of  JANNAF  Proceedings,  the  above  paper 
will  also  be  presented  at  the  AIAA  19th  Fluid  Dynamics,  Plasma  Dynamics, 
and  Lasers  Conference,  June  1987. 
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NOMENCLATURE 


A  channel  cross-sectional  area 

a  speed  of  sound 

f  frequency 

h  channel  hei ght 

M  Mach  number 

p  pressure 

R  perfect  gas  constant 

T  temperature 

t  time 

w  wave  speed 

W  wave  speed  normalized  with  average  sound  speed 

X  streamwise  coordinate;  x  ■0.  at  throat 

y  vertical  coordinate;  y  ■  0.  at  bottom  wall 
Y  ratio  of  specific  heats 

6  boundary-layer  thickness 

0  phase 

Subscri pts 
e  exit  station 

o  plenum 

t  total 

tr  transducer 

0  shock 

+  downstream- travel ing 

upstream-traveling 
*  throat 

Superscripts 

(  )'  ensemble- aver aged  pulse  fluctuation  value 
(~)  time-mean  value 
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mxiAL  Disascnxia 


16  Naval  Mx  Systens  Connaiid 
AIR-4231  (1) 

AIR-5004  (2) 

AIR-522  (1) 

AIR-524  (1) 

AIR-536  (1) 

AIR-540  (1) 

AIR-5402  (1) 

AIR-5404  (1) 

AIR-93A  (1) 

AIR-93D  (1) 

AIR-9301B  (2) 

AIR-932  (1) 

AIR-932<J  (1) 

PMA-242  (1) 

6  Naval  Sea  Systems  Coomand 
SEA-09B312  (2) 

SEA-521  (1) 

SEA-62R2  (1) 

SEA-62Z2  (2) 

1  Space  and  Naval  Warfare  Systems  Command  (SPAHAR-005) 

1  Commander  in  Chief,  U.S.  Pacific  Fleet,  Pearl  Harbor  (Code  325) 

1  Commander,  Ihird  Fleet,  San  Francisco 
1  Commander,  Seventh  Fleet,  San  Francisco 
1  Air  Test  and  Evaluation  Squadron  5,  China  Lake 
1  David  Taylor  Research  Center,  Betheada  (Code  166,  J.  Talbot) 

3  Narine  Corps  Combat  Development  Center,  Quantico  (Landing  Force  Development  Center) 
1  Naval  Air  Development  Canter,  Warminster  (Code  3014) 

1  Naval  Explosive  Ordnance  Disposal  Technology  Center,  Indian  Head 

1  Naval  Ocean  Systems  Center,  San  Diego  (Technical  Library) 

3  Naval  Ordnance  Station,  Indian  Head 

Code  5251  (1) 

Code  5253  (1) 

Technical  Library  (1) 

2  Naval  Research  Laboratory 

Laboratory  for  Computational  Physics 
Dr.  J.  Gardner  (1) 

Dr.  E.  Oran  (1) 

1  Naval  Strike  Warfare,  Fallon  (Intelligence  Library) 

3  Naval  Surface  Warfare  Center,  White  Laboratory,  Silver  Spring 

Code  K24  (1) 

Code  R16  (1) 

Code  R44,  Dr.  T.  Rsleh  (1) 

2  Naval  Technical  Intelligence  Center 

LNN,  Liaison  Officer  (1) 

MITC-60,  Library  (1) 

1  Naval  War  College,  Newport 

1  Pacific  Nlsaile  Test  Center,  Point  Nugu  (Code  1094.7,  P.  NcQuaide) 

1  Army  Amamant  Hunitlons  and  Chemical  Cnamnd,  Rock  Island  (DRSAR-LEM)  _ 

1  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground  (DRDAR-TSB-S  (STINFO)) 

2  Air  Force  Systama  Coaund,  Andrews  Air  Force  Base 

AFSC-SDZ  (1) 

AFSC-XETA  (1) 
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6  Air  Force  Armaoent  Division,  Eglln  Air  Force  Base 
AFAIL/DUG  (1) 

AFAIL/DLMA  (1) 

AFAIL/DLMI,  Aden  (1) 

AFAIL/DLO  (1) 

AFAIL/DLODL  U) 

AFAIL/FXE  (1) 

5  Air  Force  Astronautics  Laboratory,  Edwards  Air  Force  Base 
Code  NRCC  (i) 

Code  MKP  (1) 

Technical  Library  (3) 

1  Air  Force  Intelligence  Agency,  Bolling  Air  Force  Base  (AFIA/INTAW,  NAJ  R.  Esaw) 

5  Air  Force  Wright  Aeronautical  Laboratories,  Wright- Patterson  Air  Force  Base 

AFWAL/MLBX  (3) 

APVAL/POPR  (1) 

AFWAL/POPX  (1) 

2  Foreign  Technology  Division,  Wri^t-Patterson  Air  Force  Base 

FED/SDNW  (1) 

FTD/TqFA  (1) 

1  Defense  Advanced  Research  Projects  Agency,  Arlington 
12  Defense  Technical  Information  Center,  Alexandria 
1  Department  of  Defense  Explosives  Safety  Board,  Alexandria  (6-A-14S) 

1  California  Institute  of  Technology,  Pasadena,  CA  (Prof.  F.  Culick) 

1  Flow  Research  Company,  Kent,  WA  (Dr.  Wen-Huei  Jou) 

1  General  Dynamics  Corporation,  Convair  Division,  San  Di^o,  CA  (Research  Division) 

1  General  Dynamics  Corporation,  Pomona,  CA  (Via  KAVPRO) 

3  Georgia  Institute  of  Technology,  School  of  Aerospace  Engineering,  Atlanta,  GA 

Prof.  S.  Lekoudia  (1) 

Prof.  W.  Strahle  (1) 

Prof.  B.  Zinn  (1) 

1  Grumman  Aerospace  Corporation,  Betlqpage,  NY  (TIC  LOl-35) 

1  Honeywell  Incorporated,  MilitSKy  Avionics  Division,  St.  Louis  Paric,  MN 
1  Hudson  Institute,  Incorporated,  Canter  for  Naval  Analyses,  Alexandria,  VA  (Technical 
Library) 

1  Oighes  Aircraft  Co^Muiy,  Missile  Systems  Gro«;qp,  Canoga  Park,  CA  (Library) 

1  Ii^erial  College  of  Science  and  Technology,  Department  of  Mechanical  Engineering, 
London  (Prof.  J.  Whitelaw) 

3  McDonnell  Dot^las  Corporation,  St.  Louis,  MO 
Department  E2A2,  R.  Fall  (1) 

Dr.  T.  Bogar  (1) 

Dr.  M.  Sajben  (1) 

1  Raytheon  Coi^any,  Missile  Systems  Division,  Bedford,  Ma 

1  Rockwell  International  Corporation,  Missile  Systems  Division,  Duluth,  GA  (F.  Hessman) 

1  San  Diego  State  University,  San  Diego,  CA  (Prof.  N.  Nosseir) 

2  Stanford  University,  Department  of  Mechanical  Engineering,  Stanford,  CA 

Prof.  C.  Bowman  (1) 

Prof.  R.  Hanson  (1) 

1  Teledyne-Ryan  Aeronautical,  San  Diego,  CA 
1  Bie  Boeing  Company,  Seattle,  WA 

1  The  Johns  Hopkins  University,  Applied  Physics  Laboratory,  Laurel,  MD  (Kiersey) 

1  United  Technologies  Corporation,  Chemical  Systems  Division,  San  Jose,  CA  (T.  Meyers) 

1  United  Technologies  Research  Center,  East  Hartford,  CT  (Doctaaent  Control  Station,  R. 
O’Brien) 

1  University  of  California,  Berkeley,  Department  of  Mechanical  Engineering,  Berkeley,  CA 
(Prof.  J.  Daily) 

1  University  of  California,  Davis,  CA  (Prof.  C.  Law) 

1  University  of  California,  Irvine,  Department  of  Mechanical  Engineering,  Irvine,  CA 
(Prof.  W.  Sirignano) 

1  University  of  Illinois,  Department  of  Chemistry,  Chicago,  IL  (Prof.  R.  Moriarlty) 
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